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ABSTRACT 

Using reconstructed opacities, we construct solar models with low heavy- 
element abundance. Rotational mixing and enhanced diffusion of helium and 
heavy elements are used to reconcile the recently observed abundances with he- 
lioseismology. The sound speed and density of models where the relative and 
absolute diffusion coefficients for helium and heavy elements have been increased 
agree with seismically inferred values at better than the 0.005 and 0.02 fractional 
level respectively. However, the surface helium abundance of the enhanced diffu- 
sion model is too low. The low helium problem in the enhanced diffusion model 
can be solved to a great extent by rotational mixing. The surface helium and 
the convection zone depth of rotating model M04R3, which has a surface Z of 
0.0154, agree with the seismic results at the levels of 1 a and 3 o respectively. 
M04R3 is almost as good as the standard model M98. Some discrepancies be- 
tween the models constructed in accord with the new element abundances and 
seismic constraints can be solved individually, but it seems difficult to resolve 
them whole scenario. 
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Introduction 



Recent analyses of solar photospheric abundance h ave shown a significant reduction of 



the abundances of C, N, 0, and other heavy element s (ILodders 



2003 



Asplund et al.l 12004 



Allende Prieto fc Lambertl 120051 : lAsplund et all l2005t IScott et al.l 120061 ). As a result, the 
ratio of the heavy-elemen t abundance to hydrogen ab undance, Z/X, is red uced from 0.023 



(iGrevesse fc Sauval 11998. hereafter GS98) to 0.0171 flAsplund et all 12004 hereafter AGS) 
(or 0.0177 ILodders! (120031 )); and the surface Z of the Sun decreases from 0.017 to 0.0126 
(or 0.0133 ). Many solar models constructed in accord with this low Z disagree with the 
seismically inferred sou nd speed and density p rof iles, convection zone (C Z ) depth, and helium 



abundance of the CZ (IBasu fc Antia 



20041 : iBahcall et alJliooa 



Guzik et al. 



2004. 



2005|). 



Montalbdn et al. 2004; Turck-Chieze et al. 



Basu fc Antial fl2004aJ . hereafter BA04) and IBahcall et al.l (120051 ) found that the 15 - 
20 per cent increase of OPAL opacities at the base of the CZ can resolve the discrepancies 
bet ween the solar models with low Z and helioseismology. However, ISeaton &: Badnelll (120041 ) 
and iBadnell et al.l (120051 ) show ed that the inc r ease in the opacities is no more than 2.5 per 
cent near the base of the CZ. I Antia fc Basul (120051 ) and IBahcall et al.l (120051) found that 
increas i ng ne on abundance can solve the discrepancies. However. ISchmelz et al. ( 2005 ) and 
Young! (120051 ) showed that the ratio of Ne/O is consistent with the value given by AGS. 
Another attempt to resolve the problem focuses on increasing the rates of diffusion of helium 
and heavy elements. Models with enhanced diffusion rates are in better agreement with 
helioseismology than the model without the enhanced diffusion. However, the sound speed 
and the density are still far from the seismic result s, or the surface Z/X value is still too 



high and the helium abu ndance of the CZ is too low (IBasu &: Antia! l2004al iMontalbdn et al. 
20041 : iGuzik et alil2005h . 



Using helioseismic data, lAntia fc Basul (120061 ) determined t hat the solar metal abun 



dance is 0.0172 ± 0.002. This is consistent with that of GS98. iDelahaye fc Pinsonneault 



( 120061 ) found that both Fe/H and O /H are more consistent with the values of GS98 than 
the values of AGS. lAyres et al.l (120061 ) also found that the oxygen abundance is close to the 
value of GS98. If t hese resul t s are further confirmed, the problems induced by low Z will 
disappear. Howev er, lAsplundl (120061 ) argued that the AGS results are trustworthy. Recently, 
Scott et al.l (120061 ) also found that the low carbon abundance is in good agreement with the 
findings base d on e ntirely different indicators of AGS and wit h the values determ ined by 
Ireland et al.l (120061 ) from lunar grains irradiated b y solar wind . Scott et al.l (120061 ) claimed 
that their results are more reliable than those of Ayres et al. ( 20061 ). The answer to the 
problems is still an open question. 



Our aim is to construct a solar model, using recently determined abundances, which 
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can agree with seismic constraints. We apply straight multipliers to the diffusion velocity 
to enhance the rates of element diffusion. Althou gh the theoretical error of gr avitational 
settling rate is of the order of about 15 per cent (IThoul. Bahcall &: Loeblll994T). there ar e 
some significant uncertanties in the treatment of the element diffusion (IGuzik et al.ll2005l ). 
Our multipliers of the diffusion coefficients are very high, despite the fact that there is no 
obvi ous physical .justifi cation for such high multipliers, as has been pointed out by BA04 
and iGuzik et al.l (120051 ). However, the multipliers are required in our models for diffusion 
to reduce the heavy-element abundance of the CZ from the GS98 value to near the AGS 
value. In order to get the same helium in a rotating model as in a non-rotating model, a 
multiplier of the element diffusion is required because rotational mixing reduces the degree of 
gravitational sett l ing. O ther enhanced dif f usion- rate models have been discussed by BA04, 



Montalbdn et al.l (120041 ). and IGuzik et al.l (120051 ). The main difference between our models 



and those of others is that we include rotational effects in order to resolve the low heli um 
problem of the CZ. We compare our results with those of BA04 and IGuzik et al.l (120051 ) in 
Tabled 



Solar models 



2.1. Properties of our solar models 



We use the Yale Rotation Evolution Code (YREC7) to construct our solar models. 
However, we modify the code to include two tables of OPAL EOS for the Z diffusion; and we 
correct an error in rotating mo del in the calculation of the chemical compositions. We use 
reconstructed OPAL opacities (llglesias fc Rogers! Il996l ) and new low-temperature opacities 
([Ferguson et al.ll2005l ). both with the GS98 or AGS mixtures. In principle, the EOS table 
should also be reconstructed using AGS mixtures, but i t has been found that the effects of 



the change of heavy elements o n the EOS can be ignored dBasu fc Antia 



2004! . 120051 : IGuzik et al.ll2005h . We therefore use the OPAL EOS faogers et allll996f) in all 



2004a 



Bahcall et al. 



our m odels. Element diffusion is included for helium and metals ( IThoul. Bahcall fc Loeb 
1994J). Energy transfer by convection is treated according to the standard mixing- length 
theory, and the boundaries of the convection zones are determined by the Schwarzschild 
criterion. We take the solar age to be 4.57 Gyr. Luminosity L & = 3.8418 xlO 33 erg/s, and 
radius R Q = 6.9598 xlO 10 cm. 

We construct the following six models: 1) M98, a standard model with GS98 mixture 
opacities; 2) M04, a mo del with AGS mixture opacitie s; 3) M04D, same as M04 but enhancing 



the eleme nt diffusion dThou 



rotation (IPinsonneault et al 



Bahcall k Loebl ll994T ): 4) M04 R1, same a s M04D but with 



19891 ) and diffusion coefficient (Zahn 1993 ) added for shear 
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instability; 5) M04R2 and 6) M04R3, both the same as M04R1 but with different diffusion 
multipliers, as shown in Table [TJ 

Some of the parameters of the models are summarized in Table [TJ The mixing-length 
parameter a, Z init and Y init are free parameters adjusted to obtain the observed solar radius, 
luminosity and surface element abundances. Finally, R cz /R & , (Z/X) s , Y s , and Z s are the 
results of calculations at the age of 4.57 Gyr. 



2.2. Results 



Some of the calculation results are showed in Table [TJ The base of the CZ i s at 0.7335 
for M 04 and is 20 a different from the seismically inferred 0.713 ± 0.001 Rq (IBasu fc Antia 
19971 ). The surface helium abundance of 0.2294 is 6 a away from the seismically inferred 
value 0.2485 ± 0.0034 (BA04). The sound speed and density difference s between M04 an d 
the Sun are shown in Figure [TJ where those of the Sun were given by IBasu et al.l (120001 ) . 
M04 is in strong disagreement with the seismically inferred sound speed and density profiles, 
the depth of the CZ, and envelope helium abundance. 



Asplund et al.l (120041 ) suggested that increased diffusion might be able to resolve these 
disagreements. Thus, we construct the model M04D by multiplying the diffusion coefficients 
for helium and heavy elements by factors of 2.4 and 3.8 respectively. However, we have no 



physic al justification for these multipliers. This method was firstly proposed by lGuzik et al. 



( 120051 ). The b ase of the CZ of M0 4D is at 0.7168 Rq, which agrees with the seismically 
inferred value (IBasu fc Antial 119971 ) at about the 3 a level. The sound speed and density 
of M04D are better consistent with the seismic data than that of M04. The difference of 
the sound speed and density between M04D and the Sun, 5c/ 'c and Sp/p, is less than 0.005 
and 0.02 respectively. These values are close to those of the standard model M98. In fact, 
the density profile of M04D is even slightly better than that of M98. However, the surface 
helium abundance of 0.2225 is too low, which disagrees with the seismically inferred value 
0.2485 (BA04) at the level of 8 a. 



The rotational mix ing can reduce the degree of gravitational settling (jChaboyer et al. 



19951 ; lYang fc Bil 120061 ). We thus construct a rotating model, M04R1, to study the low 
helium problem of M04D. The ro tational mixing in all our models is treated as a diffusion 
process ( jPinsonneault et al.lll989l ); and the rotational mixing processes include Eddi ngton 
circu lation, Goldreich-Schubert-Fricke instability, and the secular shear instability (IZahn 
19931 ). In all our rotating models, the surface rotation rate is about 2.86xl0~ 6 rad/s at 
the age of 4.57 Gyr. The sound speed and density differences between M04R1 and the Sun 
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are shown in Figure [fl They are almost the same as those in M04D. The surface helium 
ab undance of M04R 1 is 0.2368, which is 3 a away from 0.2485 (BA04) but agrees with that 
of Kosovichev Jl997h . However, the base of the CZ is at 0.7206 R Q and disagrees with the 



seismically inferred value at about the 7 a level. 

In order to get a model which is more consistent with seismic constraints than M04R1, 
we relax the constraint of the heavy elements determined by AGS and construct models 
M04R2 with Z/X = 0.01976 and M04R3 with Z/X = 0.0208. The sound speed and density 
differences between M04R2 and the Sun are less than 0.004 and 0.01 respectively; and the 
surface helium abundance of M04R2 is 0.2454 within the constraint of observation (BA04). 
However, the depth of the CZ in M04R2 disagrees with the seismically inferred value at 
about the 7 a level. The 5c/ c and Sp/p of M04R3 are less than 0.005 and 0.02 respectively. 
The surface helium abundance and the base position of the CZ in M04R3 agree with the 
seismically inferred values at the levels of 1 a and 3 o respectively; but the envelope Z of 
0.0154 is higher than the value of AGS. This indicates that a model with high heavy elements 
is more consistent with seismic constraints than the model with AGS heavy elements. 

In the FULL1 M of BA04, the diff erences 5c/ 'c and 5p/p are larger than 0.01 and 0.05 



respectively. In the iGuzik et al.l (120051 ) models, the 5c/ c is larger than 0.005. We compare 



our results with these models in Table [TJ 



3. Discussion and conclusions 

In this study, we assume that the diffusion coefficients for helium and heavy elements 
can be enhanced respectively. Thus the discrepancies between the models with the low Z 
and helioseismological theory can be reduced, but the helium abundance of the CZ main- 
tains too low in such models. In enhanced diffusion models, there is a gradient of element 
abundances caused by element settling at the base of the C Z. The rotational secular shear 



is highly sensitive to the gradients of mean molecular weight (IPinsonneault et al.lll989l ). and 
the rotational mixing can smooth out the gradient of the element abundances. The settling 
of the helium abundance can thus be partly counteracted by the rotational mixing, and the 
surface helium abundance of the rotating model can thus be enhanced. 

After relaxing the constraint of the Z/X in model M04R3, the sound-speed profile, 
density profile, depth of the CZ, and surface helium abundance are almost the same as those 
in the standard model M98; but the surface Z of 0.0154 is higher than the value of AGS. 



The difference in the sound speed and density in our models can be improved by enhanc- 
ing the diffusion coefficients for helium and heavy elements. Thus, the maximum difference 



- 6- 



in sound speed decreases from 0.014 in M04 to 0.0045 in M04D and M04R1. The density 
profiles of the enhanced diffusion models are even slightly better than that of M98. The sur- 
face helium abundance of 0.2225 in model M04D is too low and is 8 a different from that of 
BA04. The surface h e lium abun dance of the rotatin g model M04R1 is 0.2368, which agrees 



with that of iLodderd (120031 ) and iKosovichevi (119971 ); but the posi tion of the base of t he CZ 



in this model is 7 o away from the seismically inferred position (IBasu fc Antial 119971 ). The 
surface helium abundance and the base position of the CZ in M04R3 are in agreement with 
seismic results at the levels of 1 a and 3 a respectively; but the surface heavy-element abun- 
dance of 0.0154 is higher than the value of AGS. Some problems between the new element 
abundances and helioseismological theory can be solved individually, but it seems difficult 
to resolve them whole scenario. 

We are grateful for help of Daniel Kister, useful remarks from the anonymous referee, 
and the support by the NSFC through projects 10473021 and 10433030. 



REFERENCES 

Allende Prieto, C, & Lambert, D. L. 2005, in LPICo. No. 1278, Workshop on Oxygen in 
the Earliest Solar System, 9A 

Antia, H. M., & Basu, S. 2005, ApJ, 620, L129 

Antia, H. M., & Basu, S. 2006, ApJ, 644, 1292 

Asplund, M., Grevesesse N., Sauval, A. J., Allende Prieto, C, & Kiselman, D. 2004, A&A, 
417, 751 (AGS) 

Asplund, M., Grevesesse N., Sauval, A. J., Allende Prieto, C, & Blomme, R. 2005b, A&A, 
431, 693 

Asplund, M. 2006, in IAU Symp. No. 239, Convection in Astrophysics, 9 
Ayres, T. R., Plymate, C, & Keller, C. 2006, ApJS, 165, 618 

Badnell, N. R., Bautista, M. A., Butler, K., Delahaye, F., Mendoza, C, Palmeri, P., Zeippen, 
C. J., Seaton, M. J. 2005, MNRAS, 360, 458 

Bahcall, J. N., Serenelli, A. M., & Pinsonneault, M. H. 2004, ApJ, 614,464 

Bahcall, J. N., Basu, S., Pinsonneault, M. H., & Serenelli, A. M. 2005, ApJ, 618, 1049 



-7- 



Basu, S., & Antia, H. M. 1997, MNRAS, 287, 189 

Basu, S., Pinsonneault, M. H., k Bahcall, J. N. 2000, ApJ, 529, 1084 

Basu, S., & Antia, H. M. 2004a, ApJ, 606, L85 (BA04) 

Basu, S., & Antia, H. M. 2004b, in Helio- and Asteroseismology: Towards a Golden Fu- 
ture, Proc. of the SOHO 14 / GONG 2004 Workshop, ed. D. Danesy (ESA SP-559; 
Noordwijk: ESA), 317 

Chaboyer, B., Demarque, P. & Pinsonneault, M. H. 1995, ApJ, 441, 865 

Delahaye, F., & Pinsonneault, M. H. 2006, ApJ, 649, 529 

Ferguson, J. W. et al. 2005, ApJ, 623, 585 

Grevesse, N., & Sauval, A. J. 1998, in Solar Composition and Its Evolution, ed. C. Frohlich, 
et al. (Dordrecht: Kluwer), 161 (GS98) 

Guzik, J. A., Watson, L. S. & Cox, A. N. 2005, ApJ, 627, 1049 

Iglesias, C, Rogers, F. J. 1996, ApJ, 464, 943 

Ireland, T. R., Holden, P., Norman, M. D., & Clarke, J. 2006, Nature, 440, 776 

Kosovichev, A. G. 1997, in AIP Conf. Proc. 385, Robotic Exploration Close to the Sun: 
Scientific Basis, ed. S. R. Habbal (Woodbury: AIP), 159 

Lodders, K. 2003, ApJ, 591, 1220 

Montalbdn, J., Miglio, A., Noels, A., Grevesse, N., Di Mauro, M. P. 2004, in Helio- and 
Asteroseismology: Towards a Golden Future, Proc. of the SOHO 14 / GONG 2004 
Workshop, ed. D. Danesy (ESA SP-559; Noordwijk: ESA), 574 

Pinsonneault, M. H., Kawaler, S. D., Sofia, S., & Demarqure, P., 1989, ApJ, 338, 424 

Rogers, F. J., Swenson, J., & Iglesias, C. 1996, ApJ, 456, 902 

Scott, P. C, Asplund, M., Grevesesse N., & Sauval, A. J. 2006, A&A, 456, 675 

Seaton, M. J., & Badnell, N. R. 2004, MNRAS, 354, 457 

Schmelz, J. T., Nasraoui, K., Roames, J. K., Lippner, L. A., & Garst, J. W. 2005, ApJ, 634, 
L197 



- 8- 

Thoul, A. A., Bahcall, J. N., Loeb, A. 1994, ApJ, 421, 828 
Turck-Chieze, S., Couvidat, S., Piau, L. et al. 2004, PhRvL, 93, 1102 
Yang, W. M., & Bi, S. L. 2006, A&A, 449, 1161 
Young, P. R. 2005, A&A, 444, L45 

Zahn, J. P. 1993, in Astrophysical Fluid Dynamics, Les Houches XLVII, ed. J. -P. Zahn & 
Zinn-Justin(New York: Elsevier), 561 



This preprint was prepared with the AAS IATgX macros v5.2. 



- 9 - 









1 1 1 1 1 1 I 


o 








d 




A 
















o 










d 




















O 




















X 




j 


/ / ' 




LO 






/ ■' 








\ 


/ !/? 




o 




\ ^ . 








1*7 Tt 






O 
















X 








m 






i 



I 

in 
o 



\ o 

■o 6 



in - 
o - 

?0 



0.5 




0.5 1 

Fig. 1. — A) Sound speed difference between the Sun and the model. B) Densit y difference 
betw een the Sun and the model. The solar values were derived from MDI data (IBasu et al. 
20001 ) . The long-dashed line refers to M98. The solid line indicates M04. The dash-dot-dot- 
dot line shows the results of M04D. The dotted line is given by M04R1. The dash-dot line 
and the bold solid line are corresponding to M04R2 and M04R3 respectively. 
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Table 1: Model parameters. 



Model 


1 tnit 


7 

"Wit 


a 


Multiplier 


RczJ Rq 


(z/x) s 


Y s 


Z s 


M98 


0.2794 


0.0201 


2.166 


1.0 a (1.0) b 


0.715F 


0.0247 


0.2487 


0.01809 


M04 


0.26066 


0.0148 


1.608 


1.0 (1.0) 


0.7335 


0.0174 


0.2294 


0.01322 


M04D 


0.28618 


0.019619 


1.7827 


2.4 (3.8) 


0.7168 


0.0176 


0.2225 


0.01347 


M04R1 


0.285803 


0.019497 


1.71799 


2.4 (3.8) 


0.7206 


0.0177 


0.2368 


0.01330 


M04R2 


0.2836 


0.018892 


1.67555 


2.0 (2.5) 


0.7206 


0.01976 


0.2454 


0.01462 


M04R3 


0.28358 


0.01886 


1.6889 


2.0 (2.0) 


0.7167 


0.0208 


0.2450 


0.015396 


Basu FULL1M 








1.65 


0.7233 


0.0171 


0.2244 


0.0130 


Basu FULL2M 








1.65 


0.7138 


0.0218 


0.2317 


0.01639 


Guzik model3 








3 


0.7022 


0.0196 


0.1926 


0.01552 


Guzik model5 








1.5; 4 


0.7175 


0.0206 


0.2269 


0.01561 



"The multiplier for the diffusion coefficient of the helium; 

b The multiplier for the diffusion coefficient of the heavy elements; 

c Using OPAL EOS96. liahcall et all |20o4) obtained R cz = 0.7f55 R Q . 



